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“In  an  attempt  to  understand  mobility  mechanisms,  to  predict  optimal  con¬ 
duction  conditions  for  solid  electrolytes,  and  to  develop  methods  for  calcu¬ 
lating  conductivity  in  classical  many-body  systems,  calculations  have  been 
completed  for  a  variety  of  one-,  two-,  and  three-dimensional  model  lattice 
systems.  These  calculations  involve  Monte-Carlo  simulation,  self-consistent- 
field  techniques,  and  reduced  Langevln  dynamics,  and  have  been  devoted  to  the 
actual  calculation  of  the  frequency-dependent  conductivity  tensor  oCul),  and  of 
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its  variation  with  important  experimental  variables  (concentration,  pressure, 
temperature,  trap  sites,  identity  of  mobile  ion,  lattice  geometry).  We  have 
characterized  the  relationships  among  conduction  mechanism,  carrier  density, 
ion  potentials,  and  observed  response.  We  have  examined  the  extent,  impor¬ 
tance,  and  effects  of  ion-ion  coulomb  correlations,  which  increase  conduction 
in  incommensurate  cases.  We  have  examined  the  role  of  compensators  and  the 
extent  to  which  trapping  inhibits  conduction.  The  interrelationship  of  static 
lattice  potential,  ion-ion  repulsion,  stoichiometry,  and  conduction  has  been 
articulated.  Finally,  we  have  continuously  modified  our  computational  efforts 
on  the  basis  both  of  the  emerging  mechanistic  picture  and  of  constant  inter¬ 
action  with  and  testing  by  local  experimental  groups. 

The  ambition  of  this  research  has  been  to  provide  hard  numerical  data  for 
three  purposes:  understanding  and  prediction  of  the  effects  of  variation  of 
external  parameters  (temperature,  doping,  compensation,  ion  exchange)  on  con¬ 
ductivity,  comparison  of  good  numerical  data  with  extant  formal  theory  to 
focus  on  the  most  important  variables  or  achieving  optimal  conductivity,  and 
development  of  a  mechanistic  model  (including  such  extremes  as  hopping  and 
liquid-like  diffusion)  to  correlate  the  behaviors  of  a  large  number  of  frame¬ 
work  conductors.  We  believe  that  these  ambitions  have  been  satisfied  quite 
well:  we  have  explained  extant  experiments  and  predicted  new  conductivity 
features  for  sucn\.  important  ionic  conductors  as  8”~alumina,  NaI-poly(ethylene 
oxide),  crystalline  and  glassy  0-eucryptite,  and  NASICON.  In  addition,  we 
have  developed  new  theoretical  approaches,  notably  the  trajectory-bundled 
time-dependent  self-consistent-field  method,  and  dynamic  percolation  theory, 
which  promise  to  be  very  efficient  for  planned  future  research  on  ionic  conduc¬ 
tivity  at  interfaces,  in  glasses,  and  in  alloy  solid  electrolytes. 
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Statement  of  Problem: 


Solid  electrolytes,  or  fast-ion  conductors,  are  under  intense  scrutiny  in 
a  number  of  industrial,  government,  and  academic  laboratories  for  applications 
in  fuel  cells,  in  stable,  long  shelf-life  primary  batteries,  and,  especially, 
in  high  energy-density,  stable  secondary  batteries  for  vehicular  propulsion. 
Our  research  efforts  have  been  directed  at  the  understanding  of  the  mechanisms 
by  which  conductivity  occurs  in  framework  solid  electrolytes,  and  at  possible 
means  of  enhancement  of  the  conductivity  by  selection  of  different  experimen¬ 
tal  or  performance  conditions,  involving  control  of  such  variables  as  ion  den¬ 
sity,  bulk  density,  temperature,  crystallinity,  and  counterion.  The  problem, 
then,  might  best  be  defined  as:  Mechanistic  Understanding  of  Ionic 
Conductivity,  and  Enhancement  Mechanisms  for  Ionic  Conductivity,  in  Framework 
Solid  Electrolytes. 


In  an  attempt  to  understand  mobility  mechanisms,  to  predict  optimal  con¬ 
duction  conditions  for  solid  electrolytes,  and  to  develop  methods  for  calcu¬ 
lating  conductivity  in  classical  many-body  systems,  calculations  have  been 
completed  for  a  variety  of  one-,  two-,  and  three-dimensional  model  lattice 
systems.  These  calculations  involve  Monte-Carlo  simulation,  self-consistent- 
field  techniques,  and  reduced  Langevin  dynamics,  and  have  been  devoted  to  the 
actual  calculation  of  the  frequency-dependent  conductivity  tensor  cK  m) ,  and  of 
its  variation  with  important  experimental  variables  (concentration,  pressure, 
temperature,  trap  sites,  identity  of  mobile  ion,  lattice  geometry).  We  have 
characterized  the  relationships  among  conduction  mechanism,  carrier  density, 
ion  potentials,  and  observed  response.  We  have  examined  the  extent,  impor¬ 
tance,  and  effects  of  ion-ion  coulomb  correlations,  which  increase  conduction 
in  incommensurate  cases.  We  have  examined  the  role  of  compensators  and  the 
extent  to  which  trapping  inhibits  conduction.  The  interrelationship  of  static 
lattice  potential,  ion-ion  repulsion,  stoichiometry,  and  conduction  has  been 
articulated.  Finally,  we  have  continuously  modified  our  computational  efforts 
on  Che  basis  both  of  the  emerging  mechanistic  picture  and  of  constant  inter¬ 
action  with  and  testing  by  local  experimental  groups. 

The  ambition  of  this  research  has  been  to  provide  hard  numerical  data  for 
three  purposes:  understanding  and  prediction  of  the  effects  of  variation  of 
external  parameters  (temperature,  doping,  compensation,  ion  exchange)  on  con¬ 
ductivity,  comparison  of  good  numerical  data  with  extant  formal  theory  to 
focus  on  the  most  important  variables  or  achieving  optimal  conductivity,  and 
development  of  a  mechanistic  model  (including  such  extremes  as  hopping  and 
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liquid-like  diffusion)  Co  correlate  the  behaviors  of  a  large  number  of  frame¬ 
work  conductors.  We  believe  that  these  ambitions  have  been  satisfied  quite 
well:  we  have  explained  extant  experiments  and  predicted  new  conductivity 
features  for  such  important  ionic  conductors  as  0"-alumina,  NaI-poly(ethylene 
oxide),  crystalline  and  glassy  0-eucryptite,  and  NASICON.  In  addition,  we 
have  developed  new  theoretical  approaches,  notably  the  trajectory-bundled 
time-dependent  self-consistent-field  method,  and  dynamic  percolation  theory, 
which  promise  to  be  very  efficient  for  planned  future  research  on  ionic  conduc¬ 
tivity  at  interfaces,  in  glasses,  and  in  alloy  solid  electrolytes. 
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Summary  of  Results 

We  have  obtained  significant  results  in  two  closely  related  areas: 

(a)  Mechanisms  of  Ionic  Conductivity,  in  crystalline  frameworks,  including 
possible  conductivity  enhancement  (papers  1,2,3,8,9,13). 

(b)  Development  of  Techniques  for  Calculation  of  Ionic  Conductivity  in 
non-crystalline  environments  (papers  4,5,6,7,10,11,12). 

These  efforts  have  involved  different  techniques,  personnel,  and  ambitions. 

The  mechanistic  work  (a)  largely  involved  use  of  techniques  which  we  and 
others  had  previously  developed  (notably  Langevin  dynamics  and  Monte-Carlo 
simulation)  to  study  conductivity  mechanisms  and  enhancement  both  in  such 
actual  systems  as  NASICON,  Na+  g"  alumina,  Ba++  g"  alumina  and  g-eucryptite. 

The  technique-oriented  work  (b)  was  aimed  at  extending  the  capability  for 
mechanistic  investigation  to  highly  disordered  materials,  in  particular  to  the 
glassy  framework  electrolytes  whose  development  (largely  in  British  and  French 
laboratories)  has  been  one  of  the  most  exciting  new  developments  in  the  solid 
electrolyte  area.  In  these  materials,  the  framework  lattice  is  not  periodic, 
so  that  the  treatments  of  the  framework  either  as  the  provider  of  a  periodic 
lattice  for  hopping  (on  which  the  Monte-Carlo  work  is  based)  or  as  the  source 
only  of  periodic  (static)  and  random  (dynamic)  potentials  (on  which  the 
Langevin  dynamics  work  are  based)  become  invalid.  Both  the  dynamic  percolation 
(papers  7,11)  and  self-consistent-field  (papers  4,5,6,10,12)  methods  can  treat 
this  disordered  problem,  and  we  have,  in  fact,  applied  the  former  to  polymeric 
electrolytes. 

This  summary  will  discuss  these  two  areas  of  research  independently, 
although  they  actually  comprise  two  parts  of  an  integrated  effort  on  mecha¬ 
nisms  of  ionic  conductivity  in  ordered  (a)  and  in  disordered  (b)  lattices. 
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For  convenience,  we  have  collected  the  abstracts  of  our  publications  as 
Appendix  A,  to  which  the  reader  may  refer  in  connection  with  this  summary. 

(a).  Mechanisms  of  Ionic  Conductivity,  and  of  Conductivity  Enhancement, 
in  Crystalline  Framework  Solid  Electrolytes. 

The  crystalline  framework  electrolytes  are,  technologically,  the  most 
important  class  right  now;  g  aluminas  are  the  electrolyte  for  the  advanced 
sodium/sulfur  battery.  All  of  the  crystalline  frameworks  are  characterized  by 
well-defined  conduction-path  structures,  usually  within  a  reduced  dimensionality 
(sheets  for  g-alumina,  tunnels  for  NAISCON  and  for  g-eucryptite) .  There 
exists  a  good  separation  in  time  scale  between  lattice  vibrations  (  ~10~^  sec) 
and  ionic  hopping  times  (~  10-12  sec),  so  that  Langevin  dynamics,  with  no 
memory  effects,  is  a  valid  model  for  the  ionic  motion.  The  equations  of 
motion  are  then 

mx£  *  -7i(V(xi)}  +  R(t)  -  ml**  .  (1) 

Here  the  ionic  carrier  at  position  xj  has  mass  m,  and  is  subjected  to  three 
forces:  that  derived  from  the  potential  (including  both  coulorabic  and  frame¬ 

work  terms)  and  the  random  (R)  and  damping  (  r)  forces  arising  from  the  lattice 
motion.  By  using  the  full  coulomb  raany-body  potential  in  V,  we  have  studied 
mechanism  and  dynamics  of  ionic  motion  in  tunnel  conductors  (papers  1,2, 3, 9). 

The  ionic  carriers  in  framework  solids  are  usualliy  fairly  light  (Na+,Li+) 
so  that  their  motion  is  not  overdamped.  Paper  1  uses  both  Langevin  simulation 
and  experimental  data  (absorption  spectra)  to  show  that  the  dynamics  of  these 
carriers  cannot  be  calculated  by  simply  assuming  rapid  momentum  equilibration, 
or  replacing  Langevin  by  Smoluchowski  dynamics.  It  also  defines  a  possible 
effective-potential  model  for  use  in  defining  a  simple,  one-body  picture  of 
the  carrier  dynamics  in  such  systems. 
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Paper  2  includes  an  overview  of  the  conduction  mechanisms  in  a  number  of 
solid  electrolytes.  In  particular,  it  stresses  the  substantial  differences 
expected  between  motion  in  a  commensurate  stoichiometry  (where  the  number  of 
lattice  sites  is  an  integral  multiple  of  the  number  of  carriers)  and  in  an 
incommensurate  case  (where  the  multiple  is  not  integral).  This  notion  is 
deepened  considerably  by  the  calculations  of  paper  3,  in  which  it  is  shown 
that  the  local  minima  in  the  o~p  (conductivity/density)  plot  observed  for 
commensurate  cases  always  correspond  to  hopping  behavior,  and  that  coulomb 
correlations  will,  for  this  situation,  always  reduce  the  conduction. 
Conversely,  incommensurate  systems  can,  for  strong  enough  correlations,  move 
the  ions  by  a  collective,  liquid-like  mechanism,  and  here  coulombic  repulsion 
can  actually  lead  to  increased  conductivity.  This  suggested  that  carriers 
with  high  charges  moving  in  a  framework  with  highly  incommensurate  stoichio¬ 
metry  should  conduct  well;  concomitant  experimental  work  by  Whitmore  and  his 
students  is  investigating  this  situation. 

Paper  9  applies  the  Langevin  technique  to  study  carrier  trapping  by 
counterions  or  by  crystal  defects.  We  show  that  the  nonperiodic,  deep-well 
potentials  due  to  these  entities  compete  with  the  periodic  framework  potential 
and  with  the  interionic  coulomb  repulsion  to  produce  very  unexpected  behavior. 
By  trapping  one  or  more  carriers,  these  deep-well  potentials  can  render  a  for¬ 
merly  commensurate  density  effectively  incommensurate,  thus  actually 
increasing  the  overall  conductivity  by  trapping  some  carriers!  For  incommen¬ 
surate,  good  conductor  situations,  on  the  other  hand,  the  trapping  does  indeed 
reduce  the  conductivity. 

The  Langevin  dynamics  method  has  several  advantages;  in  particular,  it 
gives  the  conductivity  for  all  frequencies,  rather  than  just  at  DC.  Our 
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studies  have  thus  been  useful  for  comparison  with  optical  and  microwave 
spectra,  as  well  as  with  overall  dc  conduction,  and,  as  such,  have  helped 
clarify  such  concepts  as  attempt  frequencies,  correlation  factors,  and  collec¬ 
tive  effects. 

For  the  particular  case  of  g"  aluminas,  there  has  been  a  great  deal  of 
experimental  study  of  conductivity  under  a  variety  of  conditions.  For  this 
particular  system,  moreover,  the  hopping  model  is  well  justified  by  the  height 
of  the  barrier  and  the  vacancy  mechanism  of  dc  transport.  Accordingly,  we 
have  used  a  generalized  form  of  the  Monte-Carlo  method  to  analyze  conductivity 
in  several  g”  alumina  systems.  Paper  8  considers  the  quasiparticle  descrip¬ 
tion  of  conduction  near  the  half-filled  lattice  case  (p  ~0.50)  appropriate 
for  Pb++  or  Ba++  or  Ca++  g"  aluminas.  We  find  two  "superdefects"  which 
dominate  the  conduction.  Both  form  on  the  superlattice  which  consists  of  the 
alternate  sites  (every  second  site)  of  the  physical  g"  alumina  Beevers-Ross 
lattice.  At  low  temperatures,  the  "superinterstitials"  dominate  conduction  for 
p  >  0.50,  "supervacancies"  for  p  <  0.50.  While  the  latter  exhibit  low  mobility, 
the  former  are  very  mobile,  and,  accordingly,  the  optimal  stoichiometry  for 
highest  conductivity  increases  from  p  =  0.50  at  very  high  temperatures  to  p  = 
0.67  at  low  temperatures,  as  correlation  effects  become  important  and  super¬ 
interstitial  hops,  rather  than  independent  ion  hops,  provide  the  effective 
step  for  conductivity.  In  further  work  to  be  published,  we  extend  this  notion 
to  the  analysis  of  low-angle  scattering  data,  ordering  phase  transitions,  and 
relaxation-time  behavior  in  these  g"  aluminas. 

Several  experimental  groups  have  prepared  g"  alumina  structures  with  two 
different  mobile  species  present  in  the  conduction  plane.  In  paper  13,  we 
employ  the  Monte-Carlo  technique  to  analyze  the  expected  conductivity  phase 
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diagram  for  the  Ba++/Na+  g"  alumina  case.  We  find  that,  as  the  Ba++  ion 
exchange  process  begins,  Na+  conduction  first  increases  as  the  number  of 
available  vacancies  grows,  then  decreases  due  to  strong  correlation  effects. 
Once  again,  a  point  of  optimum  conductivity  is  found,  and  once  again  it  is 
temperature-dependent.  For  reasonable  temperature  (T  =  300°C),  the  conduc¬ 
tivity  maximizes  for  an  ion  exchange  of  about  30%  of  the  initial  Na+  ions 
present;  again,  experimental  work  is  under  way  in  Whitmore's  group  to  investi¬ 
gate  this  prediction. 

These  mechanistic  studies  have  accomplished  three  aims:  they  have  demon¬ 
strated  that  Langevin  dynamics  and  hopping  models  are  useful  techniques  for 
studying  ionic  conductivity  (and  that  Smoluchowski  dynamics  is  less  useful), 
they  have  furnished  mechanistic  understanding  of  the  published  experimental 
work  on  framework  electrolytes,  and  they  have  predicted  conditions  under  which 
maximal  conductivity  can  be  achieved  for  these  materials. 

(b)  Development  of  Techniques  for  Calculations  of  Ionic  Conductivity  in 
Non-Crystalline  Environments 

The  crystalline  framework  conductors  described  above  have  dominated 
experimental  work  on  fast-ion  conductors  until  quite  recently.  Since  1981, 
however,  there  has  been  growing  interest  and  activity  in  the  area  of  polymeric 
and  glassy  solid  electrolytes.  In  these  materials,  both  the  mobile-ion  array 
and  the  underlying  framework  lattice  are  highly  disordered.  For  such  materials 
the  Langevin  and  hopping  models  become  difficult  to  apply.  Accordingly,  we 
have  been  developing  new  theoretical  methods  to  deal  with  these  systems. 


The  essential  nature  of  ionic  motion  in  such  glassy  framework  conductors 
as  glassy  g-eucryptite  or  NASIGLAS  or  Li20/Br2C>3/Si02  involves  coulombically- 


correlated  mobile  ions  moving  in  a  highly  structured  but  disordered  back¬ 
ground.  The  theoretical  approach  for  study  of  such  a  system  would  involve 
first  accounting  for  the  coulombic  correlations  to  define  an  effective  one- 
particle  static  or  dynamic  potential  for  the  ionic  motion,  and  then  to  embed 
the  resulting  single  effective-particle  dynamics  in  the  disordered  potential 
of  the  glassy  framework.  We  have  taken  a  significant  step  toward  the  develop¬ 
ment  of  such  a  one-particle  effective  potential  by  means  of  the  static  and 
time-dependent  semiclassical  self-consistent-field  (SCF)  method,  which  is  pre¬ 
sented  and  applied  to  some  small-system  test  cases  in  papers  4,5,6,  and  12. 

The  idea  of  these  SCF  studies  is  a  simple  one.  If  the  wavefunction  for  a 
system  of  N  particles  is  approximated  as  a  simple  product: 

f(x1,x2...xn)  =  n  <H  (xO 


with  4>  the  single-particle  functions  and  Xi  the  coordinates  of  particle  i, 
then  the  variational  principle  yields  Hartree-like  equations  of  motion  for  the 
individual  particles 

[hi(xi)  -  ei(xi) ]  fa.(xi)  =  0  ,  (3) 

where  is  the  energy  of  the  ictl  particle,  and  the  effective  hamiltonian  is 
given  by 


hi  =  Tt  +  <  IT  <^j  (xj  )|  VTOt|  IT*  $j(xj)  >  , 


j 


(4) 


with  T^  the  kinetic  energy  of  the  ic^  particle.  Thus  each  particle  moves  in  a 
field  defined  by  the  average  of  the  other  particles.  If  ^  is  taken  as  time 
dependent,  then  (3)  becomes  a  dynamical  result: 

ift  3ij>i(xi,t) 


at- 


hi(xi,t) ^(xi.t) 


(5) 
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In  papers  4,  6  and  12  we  apply  the  static  SCF  result  (4)  to  the  analysis 
of  vibrational  states  for  three  coupled  particles.  We  develop  the  semi- 
classical  techniques  needed  to  solve  (4)  for  realistic  systems,  and  consider 
the  choice  of  coordinate  systems  in  which  the  problem  is  best  solved.  Paper  6 
is  of  particular  interest  for  applications  to  framework  glasses,  since  it 
analyzes  the  occurrence  of  structural  instabilities  and  ordering  due  to  vibra¬ 
tional  motion;  similar  ordering  phenomena  are  observed  in  frameworks,  and  it 
will  be  of  real  interest  to  see  if  they  survive  in  the  disordered  background 
potential  of  a  glass.  This  paper  also  deals  with  the  possible  definition  of 
collective  coordinates  for  problems  involving  a  large  number  of  modes;  analy¬ 
sis  in  terms  of  such  coordinates  might  prove  particularly  useful  in  extended 
systems  such  as  those  in  glassy  electrolytes. 

In  some  disordered  electrolytes,  conductivity  beomes  efficient  only  above 
the  glass  transition  temperature;  these  systems  include  polymeric  electrolytes 
based  on  poly(ethylene  oxide)  salt  complexes.  For  these  materials,  the  entire 
idea  of  hopping  on  a  quasistatic  lattice  is  invalid,  since  the  disorder  of 
the  host  material  in  which  the  ions  move  is  dynamic  rather  than  static.  Then 
the  ordinary  theoretical  tools  used  to  study  ionic  motion  will  all  fail,  since 
the  entire  lattice  becomes  dynamically,  rather  than  statically,  defined.  The 
only  approach  to  this  problem  which  has  been  used  in  the  literature  is  based 
on  quasithermodynamic  constructs,  notably  the  free-volume  and  configurational- 
entropy  models.  We  have  developed  a  microscopic  picture  for  ionic  motion  in 
these  materials,  based  upon  a  dynamic  percolation  model  for  the  ionic  motion. 
The  model  is  characterized  by  two  timescales,  one  defining  the  time  for  the 
ion  to  jump  from  one  coordination  geometry  to  another,  and  the  second  corre¬ 
sponding  to  the  time  for  configurational  changes  to  occur  in  the  polymeric 
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framework;  this  latter  is  called  the  renewal  time,  and  glasses  and  ceramics 
are  typified  by  an  infinite  renewal  time. 

In  paper  7,  we  compare  the  predictions  of  dynamic  percolation  theory  with 
free-volume  ideas;  the  free  volume  is  in  fact  correlated  with  the  renewal  time, 
and  zero  free  volume  corresponds  to  curtailed  configurational  freedom,  to 
infinite  renewal  time,  and  to  zero  conductivity.  The  dynamic  percolation 
picture  is  richer  and  more  reflective  of  the  real  systems,  however,  since 
there  are,  in  fact,  two  microscopic  motion  parameters  in  these  systems.  In 
particular,  the  kinetic  effects  of  free-volume  motion  are  present  in  the 
percolation,  but  not  in  the  free-volume,  formulations.  This  becomes  important 
when  long-range  motions  such  as  dc  conductivity  are  considered,  since  the 
moving  ion  must  see  the  open  channel  in  the  direction  in  which  it  wants  to 
move  rather  than  in  any  arbitrary  location  in  the  lattice.  This  notion  of 
correlated  lattice  motion,  a  natural  one  in  the  dynamic  percolation  scheme,  is 
entirely  absent  in  the  (statistical)  free-volume  or  configurational  entropy 
pictures. 

Finally,  we  have  examined  the  frequency-dependence  of  the  ionic  conduc¬ 
tivity  in  the  dynamic  percolation  model  (paper  11).  This  is  of  real  interest, 
since  the  nmr,  optical,  and  microwave  investigations  of  these  materials  have 
been  among  the  most  complete  experimental  studies,  and  since  the  quasithermo¬ 
dynamic  approaches  are  inherently  unable  to  deal  with  any  dynamic  phenomena. 

We  find  several  interesting  formal  properties  of  the  dynamic  percolation 
picture;  one  of  the  most  intriguing  is  a  convincing  proof  that  the  renewal 
frequency  y  and  the  applied  field  frequency  u  enter  into  the  final  conduc¬ 
tivity  expression  only  in  the  combination  (o>  +  iy).  Thus  renewal  of  the 
lattice  affects  the  hopping  ions  like  an  applied  field  of  complex  frequency. 


We  can  also  show  that  even  in  reduced  dimensionalities  in  which  static  perco¬ 
lation  is  a  finite-length  process  (that  is,  no  diffusion  occurs  below  the 
percolation  threshold),  the  dynamic  percolation  picture,  in  agreement  with 
experiment,  always  exhibits  diffusive  behavior,  with  the  mean-square  displace¬ 
ment  a  linear  function  of  time.  This  dynamic  percolation  model  is,  so  far  as 
we  are  aware,  the  only  microscopic  technique  for  studying  ionic  conductivity 
in  dynamically-disordered  materials,  and  we  intend  to  pursue  further  its 
implications  for  mechanism,  for  optical  conductivity,  and  for  the  importance 
of  local  trap  sites. 


These  studies  of  methods  for  dealing  with  highly-disordered  electrolytes 
should,  we  feel,  be  as  useful  and  helpful  for  dealing  with  some  of  the  very 
newest  materials  (such  as  polymeric  and  glassy  electrolytes)  as  the  Langevin 
and  Monte-Carlo  methods  have  proved  in  study  of  the  crystalline  framework, 
materials.  We  thus  feel  that  this  project  has  provided  both  numerical  and 
mechanistic  insights  and  some  methodological  advances  for  the  study  of  solid 
electrolytes. 
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V.. rational  Levels  In  the  Self-Consistent-Fleld  Approximation  with  Local  and  Normal 
Modes.  Results  for  Water  and  Carbon  Dioxide 
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Calculations  are  presented  for  vibrational  energy  levels  of  nonbending  HjO  (at  correct  and  linear  geometries) 
and  for  linear  HDO,  DTO,  0**04,  and  C**0**0.  Results  obtained  from  accurate  numerical  studies  are  compared 
to  those  generated  by  using  the  self-consistent-field  (SCF)  approximation  in  both  normal  and  local  coordinates. 
We  find  that  for  all  cases  except  H20,  SCF  corrections  systematically  and  considerably  improve  the  zero-order 
uncoupled  normal-mode  eigenvalues  and  render  them  superior  to  the  local  modes  even  for  part  of  the  vibrational 
overtone  region  of  HDO  and  DTO.  For  H20,  the  local-mode  picture  remains  better  than  SCF,  due  partly  to 
the  small  size  of  the  Wilson  coupling  in  local  modes  (or  equivalently  to  the  near-degeneracy  of  the  normal  modes). 
These  results  demonstrate  the  quality  of  the  SCF  as  the  best  independent-mode  approximation.  Moreover, 
unlike  the  cruder  deooupled-mode  results,  the  SCF  frequencies  are  in  many  ernes  of  nearly  spectroscopic  accuracy. 


I.  Introduction 

The  use  of  a  local-mode  description  of  vibrational  states 
was  developed  by  Henry  and  Siebrand,  first  in  the  context 
of  radiationless  decay1  and  later  for  characterization  of  the 
overtone  region  of  X— H  (X  =  0,  C,  N)  stretches  in  several 
molecules.2  More  recently,  a  number  of  workers  have  used 
local  modes  for  both  theoretical  and  experimental  study 
of  the  X---H  overtone  region;  we  cite  specifically  the  ele¬ 
gant  work  of  Swofford,  Long,  and  Albrecht; 3  extensive 
applications  and  analysis  have  been  presented  by  Sage44 
and  by  Wallace.*  Although  the  local-mode  picture  is 
principally  used  in  the  overtone  region,  Gelbart  and  co- 
workers  have  pointed  out  that,  in  certain  instances,  even 
the  fundamental  absorption  line  may  be  best  characterized 
by  local  modes.*  More  formal  studies  of  local  modes  and 
their  relation  to  normal  modes  were  carried  out  by  Bnimer 
and  Jaffa,7  who  showed  that,  in  a  classical  two-mode  model. 
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a  rigorous  distinction  can  be  drawn  between  local-  and 
normal-mode  behavior  based  on  the  characteristics  of 
Poincare  surfaces  of  section.  Lawton  and  Child/  using 
classical  trajectories,  have  studied  the  local-  and  normal- 
mode  behavior  of  H20  and  observed  the  transition  from 
normal  to  local  as  energy  increases. 
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F.  Williams,  ibid.,  49.  1860  (1968). 
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Dissociation  dynamics  of  vibratlonaily  excited  van  der 
Waals  clusters:  l2XY  l2+X+Y  (X,  Y  =  He,  Ne) 
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The  dynamics  of  sequential  dissociation  processes  of  the  type  Xl:(v)Y  —  X  +  l.(i>')Y  — »X  +  Y  +  MO 
(X,Y  —  Ne,  He)  are  studied  usie|  classical  trajectory  calculations  and  a  recently  presented  classical  version  of 
the  time-dependent  self-consistent  field  (TDSCF)  method.  The  results  obtained  indicate  the  presence  of 
significant  dynamical  correlation  effects  of  the  rare-gas  atoms  on  each  other  despite  the  negligible  direct 
interaction  between  them;  this  is  in  qualitative  agreement  with  experimental  findings.  Good  agreement  is 
found  for  the  rate  constants  and  the  variation  with  rare  gas  as  well  as  the  branching  ratio* 

(NeljHe— *Nel,  +  He  vs  — *  I, He  +  Ne)  calculated  from  TDSCF  and  from  classical  trajectories.  Both 
classical  trajectories  and  TDSCF  show  an  essentially  impulsive  dissociation  mechanism,  in  which  dissociation 
typically  follows  a  considerable  aumber  of  vibrations,  and  is  due  to  a  relatively  rare  internal  hard  collision 
between  an  I  atom  and  the  rare  gas.  As  in  the  three-body  I.X  case,  this  mechanism  differs  from  that  in  the 
RRKM  strong  coupling  model  Energy-  and  momentum-gap  relations,  based  on  the  weak-coupling  picture, 
are  found  to  be  relatively  successful  but  fail  to  describe  the  dynamics  quantitatively. 
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I.  INTRODUCTION 

Among  the  most  detailed  measurements  of  the  rates 
of  unimolecular  decay  of  molecules  are  the  classic  ex¬ 
periments  of  Levy  and  co-workers1-5  on  the  vibrational 
predissociation  of  van  der  Waals  clusters  involving  rare 
gases  and  I2.  For  Helt,  for  example,  the  experiment 
involves  the  following  three  steps: 

HeI1(X12)  +  /iir-HeIa(i>\  B*^)  .  (1) 

Hel,(t/,  S)-He+II(«/'f  fl) ,  (2) 

M*',  fi)-I,(X)  +  A*'  .  (3) 

By  measuring  the  absorption  llnewidth  associated  with 
Eq.  (1),  the  rate  of  decay  of  the  initially  prepared  com¬ 
plex  is  determined.2  Measurement  of  fluorescence  in 
step  (3)  provides  vibration -rotation  distributions  asso¬ 
ciated  with  decay  of  Hel,. 1  Since  the  initial  state  of  tf 
may  be  varied  simply  by  changing  the  excitation  fre¬ 
quency  v,  it  is  possible  to  map  out  both  the  decay  rate 
and  It  vibrational  distribution  as  a  function  of  initial 
energy. 

In  addition  to  the  triatomic  complexes  Helx,  Nelx, 

Arlx,  etc.,  a  number  of  tetratomic  and  larger  complexes 
have  been  studied3-*,  including  He,It,  Ne2I2,  HeNel,, 
ArHelt,  and  Artlt.  In  this  case,  the  absorption  line- 
width  measures  the  rate  of  the  first  atom  dissociated, 
and  the  Ix  fluorescence  determines  the  It  vibrational 
distribution  after  all  rare  gas  atoms  have  dissociated. 
Measurements  on  tetraatomic  clusters  have  made  it 
possible  to  study  such  questions  as  the  Importance  of 
correlation  between  different  dissociative  steps,  and 
the  nature  at  energy  migration  between  different  bonds. 
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The  wealth  of  data  provided  by  these  experiments  has 
prompted  several  theoretical  studies  of  vibrational  pre¬ 
dissociation,  all  of  which  have  thus  far  been  directed  at 
triatomic  dissociation  using  reduced  dimensionlity 
models  (either  collinear  or  T-shaped).  These  include 
the  early  work  of  Child,  *  the  coupled  channel  and  dis¬ 
torted  wave  studies  of  Beswick  and  Jortner, 3  the  Golden- 
Rule  model  of  Ewing,  *  the  quasiclassical  trajectory  cal¬ 
culations  of  Woodruff  and  Thompson,  *'*'  the  perturbed- 
wave-packet  study  of  Vlswanathan  et  al. ,  *<w  and  the 
time-dependent  self -consistent-field  treatment  of  G<v  oer 
et  al. 10  While  some  of  these  theories  are  not  intend  - 
to  be  quantitative,  it  is  significant  to  note  that  there  -  * 
good  agreement  between  the  lifetimes  predicted  by  the 
coupled  channel  and  the  trajectory  methods  for  colllr.»ar 
Hel*.  *  Agreement  of  the  Helz  coupled  channel  rate  r  - 
stants  with  experiment  is  also  pretty  reasonable7"”  . 
the  T-shaped  model  and  a  potential  composed  of  sums 
of  Morse  functions.  Of  course  it  should  be  remembered 
that  the  parameters  in  the  potential  are  not  known  inde¬ 
pendently,  and  thus  are  adjusted  to  fit  experiment  in 
making  comparisons.  Other  theoretical  studies  of 
vibrational  predissociation  in  van  der  Waals  molecules 
include  Cl*-Cl2  by  Brady  ef  nl. 11  and  Ar-HCl  by  Chiidu 
and  Dunker  and  Gordon, 13  and  Ar  •••  Nj,  Helt  by 
Shapiro. 14  In  addition,  Rice  and  co-workers**  have 
studied  He  + 12  collisions  at  energies  only  slightly  above 
the  dissociation  threshold.  In  none  of  this  work  has  the 
competition  and  correlation  between  different  decay 
channels  in  a  polyatomic  dissociation  process  been 
examined. 

In  thts  paper,  we  use  a  combination  of  quasiclassical 
trajectory  methods  and  time -dependent  self -consistent- 
field  (TDSCF)  methods1*  to  study  vibrational  predisso¬ 
ciation  in  collinear  models  of  the  tetra-atomic  clusters 
HexIt,  Ne,lt,  and  HeNelj  and  the  triatomic  clusters  Helg 
and  Neff.  The  Isomers  He-It-Ne  and  HeNe-It  are  both 
considered.  There  are  two  purposes  for  our  calcula¬ 
tions:  First,  we  are  concerned  with  mechanistic  ques¬ 
tions,  such  as  the  applicability  of  simple  impulsive 
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Structural  instabilities  in  molecular  clusters  have  been  observed  in 
several  recent  classical  simulations.  They  arise  from  changes  in  the  dyna¬ 
mical  interactions  among  the  vibrational  modes  as  energy  increases  and  are 
entirely  unrelated  to  vibration  ally-induced  instabilities  of  Jahn-Tcller  type. 
We  show  that  the  self-consistent  field  (SCF)  approximation  for  vibrations 
provides  a  simple  predictive  model  for  such  instabilities,  in  which  the  onset 
of  instability  is  shown  to  arise  from  changes  in  the  effective  (mean)  potential 
acting  on  each  mode,  the  total  energy  changes.  Exact  classical  dynamics 
was  used  to  test  the  SCF  prediction  for  a  simple  two-mode  model,  with  the 
result  that  the  instability  onset  is  correctly  predicted  by  the  SCF  approxi¬ 
mation.  We  comment  on  the  relevance  of  these  mean  potentials  for  charac¬ 
terizing  related  phase  transitions  in  dusters. 


1.  Introduction 

There  is  considerable  current  interest  in  the  question  of  structural  in¬ 
stabilities  in  molecular  cluster  systems,  with  special  attention  focused  on  the 
relationship  of  such  instabilities  to  phase  transitions  in  the  macroscopic  limit. 
Nearly  all  the  evidence  for  such  transitions  comes  from  classical  numerical 
simulations,  using  both  Monte  Carlo  and  molecular  dynamics  methods,  on 
molecular  clusters.  For  example,  Etter  et  al.  [1]  have  investigated  (C02)„ 
for  2$n^l3.  They  observed  that  even  for  «=* 2,  the  cluster  shows  an  in¬ 
stability  (or  structural  change)  as  a  function  of  energy  (or  of  temperature)  ; 
the  instability  becomes  more  pronounced  as  the  size  n  increases.  A  qualitatively 
similar  pattern  was  observed  in  the  work  of  Nauchitcl  and  Pertsin  [2]  on  noble- 
gas  clusters.  They  commented  on  the  possible  relationship  of  such  structural 
changes  to  solid-solid  and  solid-liquid  phase  transitions  in  the  corresponding 
macroscopic  materials.  An  interesting  paper  by  Rao,  Berne  and  Kalos  [3] 
studied  both  thermodynamic  and  structural  properties  of  finite  clusters  in  a 
sea-ch  for  transition  behaviour  ;  their  figures  show  striking  changes  of  orienta¬ 
tional  arrangement  as  a  function  of  energy.  Earlier  investigations  [4]  focused 
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Polymeric  solid  electrolyte*  offer  a  difficult  problem  from  the  viewpoint  of  under* 
standing  the  charge  transport  mechanism.  While  quaaltkermodynamle  theories  (config¬ 
urational  entropy,  free  volume)  are  useful  for  rationalizing  the  behavior  of  these 
materials,  they  do  not  really  amount  to  a  microscopic  picture.  We  have  developed  a 
dynamic  bond  percolation  (DBP)  model  to  describe  Ionic  conductivity  In  these 
materials.  The  DBP  model  Is  based  on  a  master  equation  describing  ton  hop*  among 
sites.  The  percolation  aspects  are  Included  by  making  the  bonds  between  site*  ran¬ 
domly  open  or  closed.  The  dynamical  aspect  Is  due  to  the  configurational  motions  of 
the  polymer,  and  results  In  the  variation  of  the  bond  assignments  as  open  or  closed. 
The  relationship  of  DBP  to  free  volume  theory  Is  sketched;  this  Involves  a  specific 
consideration  of  kinetic  effects  on  free-voluaa  motion. 


I.  Introduction 

Since  their  discovery  by  Wright  ( 1 |  and  early 
exploration  by  Aroand  [2],  polymeric  solid 
electrolytes  have  been  the  focus  of  Intense 
Investigation  In  several  laboratories.  Many 
Interesting  observations  have  been  reported. 
Including  vlbraelonal  spectra  leading  to  a 
solvation  scheme  for  the  eatloaa  [3),  thermo- 
chemical  studies  implying  a  salt  lattice  energy 
threshold  for  complex  formation  12,6],  network 
studies  shoving  the  Increase  of  conductivity 
with  decreasing  glass  transition  temperature 
1 5 1 ,  electrochemical  studies  showing  both 
negligible  electronic  conductivity  (6]  and 
significant  anionic  and  cationic  motion  (7|, 
microwave  conductivity  data  Indicating  simi¬ 
larity  of  free  polymer  and  polymer-salt  complex 
response  for  freauendes  above  0.5  GHz  (8!, 
morphological  and  thermal  analysis  studies 
Indicating  mesophase  behavior  19],  activation 
studies  showing  an  Important  rola  for  salt-frea 
polymer  In  the  conductivity  [10|,  viscoelastic 
studies  showing  an  lncrsase  in  glass  transition 
temperatur*  with  salt  concentration  [2, 11 1,  and 
a  number  of  new  polymer  hosts  (11,12!  In  addi¬ 
tion  to  the  original  materials,  polyethylene 
oxide  (PEO)  and  polypropylene  oxide.  The 
Interpretation  of  the  conductivity  In  these 
largely  amorphous  polymarlc  materials,  which 
are  generally  studied  well  above  their  glass 
transition  temperatures  Tg,  la  nearly  always 
given,  following  Arnand  [13],  la  tarns  of  the 
Vogel-Tammann-Fulcher  (VTF)  equation  of  glass 
science  ! 14]  or,  more  extensively,  in  terms  of 
tha  free  volume  theory  (FVT)  of  Cohen  and 
Turnbull  (15,16). 

The  VTF  equation  ltaelf  Is  of  smplrlcsl  origin, 
though  It  esn  b*  derived  (IS)  straightforwardly 
In  the  context  of  FVT.  The  VTF  equation  for 
th*  conductivity  o  of  a  polymer/sslt  complsx 
may  bo  wrlttsn 


a  -  AoT**/2  «*EA/k(T~T0 >  (I) 

Hart  Ao  is  a  (weakly)  temperature  dependent 
prefactor,  and  To  Is  the  so-called  equilibrium 
glass  transition  temperatura  (17,18|;  Eg  has 
tha  form  of  an  activation  energy.  Equation  (1) 
pradlcts  precisely  the  type  of  curved  activa¬ 
tion  plot  (In  oT  vs  T“l )  actually  observed  for 
many  [13]  of  thosa  complexes,  though  hysteresis 
[111  and  complications  dua  to  uncomplexed 
polymer  [10)  may  be  observed  In  certain  temper¬ 
ature  and  stoichiometric  regimes.  Bscaus* 
eq.  (1)  may  be  derived  from  FVT  [151,  and  slnca 
It  characterlzas  tha  conductivity  of  most  PEO- 
salt  complexes  fairly  satisfactorily,  it  la  now 
quite  standard  to  use  FVT  to  discuss  tha  Ionic 
conductivity  of  polymeric  solid  electroltytas. 
Thera  are,  however,  several  drawbacks  to  the 
FVT  approach.  Most  importantly.  It  la  not 
Intended  to  apply  to  the  diffusion  of  small 
Iona  within  a  solid  comprised  of  much  largar 
polymeric  chain  segment*  (nor,  originally,  to 
polymers  themselves,  for  which  a  configura¬ 
tional  entropy  model  [19-22)  Is  In  several 
senses  preferable);  rather.  It  Is  used  to 
describe  viscosity,  thcrmodyanalc  data,  and 
relaxation  behavior  of  glass-forming  materials 
[15,16,23-25).  While  Its  extension  to  neat 
polymars  may,  with  caution,  be  carried  out.  Its 
extension  to  lonle  motion  within  polymars  la 
more  problematic.  In  particular,  one  of  the 
fundamental  assumptions  of  FVT,  th*  fadla  and 
rapid  redistribution  of  the  free  volume,  may 
fall  to  hold  In  the  case  of  polymer  electro¬ 
lytes  with  low  fluidity. 

We  have  recently  presented  a  microscopic  theo¬ 
retical  treatment  of  Ionic  motion  In  polymere 
[26|,  baaed  on  a  dynamic  bond  percolation 
theory  (DBP)  model.  This  nodal  take*  cogni¬ 
zance  of  th*  fact  that  Ionic  hopping  can  occur 
on  time  scales  faster  than  typical  reorganlxa- 
tlon  time*  of  ehe  polymeric  host,  and  that  tha 
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Superdefects  In  Na  0*  Alualna:  Coapucer  Simulation  of  Ionic  Conductivity  and 

Conduction  Mechanism* 
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Mobtle-lon  correlation  phenomena  In  Na+~S  "  alualna  are  studied  using  the  Monte-Carlo 
(MC)  method  to  analyze  a  lattice  hopping  model  of  ionic  behavior.  Ionic  conductivity 
is  calculated  for  various  lon-slte  occupancies,  p's,  and  a  wide  range  of  temperaturee. 
Significant  curvature  of  the  Arrhenius  plot,  which  resembles  certain  experimental  data 
on  s”  aluminas.  Is  found  for  p  *  0.5.  This  observation  Is  analyzed  In  terms  of  for¬ 
mation  of  a  long-range  ordered  superlattice.  The  mechanism  of  charge  transport  on  the 
superlattice  Involves  movement  of  eoaplex  defects  of  the  vacancy-type,  which  we  call 
supervacancies,  or  of  the  lnterstltlal-type,  which  we  call  superlnterstltlals.  The 
Nernst-Elnsteln  relation  holds  for  superdefects.  Finally,  the  opdaum  Ion  density  for 
conduction  Is  estimated  for  different  temperatures. 


I.  Introduction 

One  particularly  Important  and  difficult  Issue 
In  understanding  mechanisms  for  motion  of  the 
mobile  Ions  In  framework  Ion  conductors 
Involves  the  role  of  lon-lon  correlations  In 
the  conduction  process  (1-6).  Comes  et  al.  IS] 
have  shown  from  X-ray  scattering  Investigations 
that  ordering  of  the  Na*  Ion*  and  vacancies  In 
che  conduction  slab  of  Na  8’  alumina  sets  la  at 
lower  temperatures  for  composition*  In  which 
che  site  fraction,  p,  occupied  by  Na*  is  dose 
to  0.50.  Such  studies  have  raised  a  number  of 
question!  about  conduction  mechanism*  In  such 
structures.  (1?  Over  what  range  of  p  can 
ordering  effects  be  observed?  (2)  What  Is  the 
dependence  of  Ionic  conductivity  on  composition 
and  on  temperature  In  and  near  this  range? 

(3)  What  Is  the  optimal  stoichiometry  for  maxi¬ 
mum  ionic  conductivity?  (A)  What  1*  the 
“a*  ion  conduction  mechanism  In  such  structure* 
and  how  does  the  mechanism  change  a*  a  function 
of  p? 

We  have  attempted  to  answer  these  questions  by 
means  of  a  computer  simulation  study  of  Na*  ion 
conduction  In  S"-alumlna  as  a  function  of  tem¬ 
perature  and  crystal  composition  using  the 
well-known  Monte  Carlo  (MC)  method  (6-9J. 
Although  this  method  suffers  some  disadvantages 
when  compared  with  the  molecular  dynamics 
method  for  studying  frequency-dependent  conduc¬ 
tion  (10, 11 1,  nevertheless  It  hes  proved  quite 
valuable  In  developing  a  better  understanding 
of  the  d.e.  conduction  behavior  In  this  class 
of  compounds .  Indeed,  the  pioneering  MC  studies 
of  Murch  and  Thorn  (7]  on  Ion  conduction  In 
8  and  8“  aluminas  have  been  especially  helpful 
In  elucidating  the  role  of  correlation  effects 
on  the  conductivity  of  these  materials.  Other 
anelytlc  schemes  have  also  beam  devised  for  the 
Investigation  of  Ionic  transport  In  framswork 
conductors.  In  particular.  Sate  and  Klkuchi 
112-14)  have  developed  a  path-probability 
method  (PPM)  which  Is  a  type  of  duster- 


variation  method  based  on  a  kinetic  fora  of  the 
quael-chealeal  approximation.  Some  discrepan¬ 
cies  have  arisen  between  the  predictions  of  the 
MC  method  and  those  of  the  PPM  for  NaB*  alumina, 
especially  In  the  compositional  region  near  p  * 
0.50;  one  goal  of  the  present  study  la  to 
clarify  this  discrepancy. 

Perhaps  the  most  Interesting  result  of  our  study 
Is  the  prediction  of  what  we  term  'super- 
defects*  (supervacancle*  and  superlnter¬ 
stltlals).  At  lower  temperatures  an  ordered 
superlattice  exists  at  the  exact  composition 
p  "  0.50;  then  small  deviations  from  the  ideal 
composition  p  -  0.50  result  In  che  Introduction 
of  Interstitials  or  vacancies,  which  arrange 
themselves  Into  "superdefects*  on  che  euper lat¬ 
tice.  Such  superdefects  move  quite  readily 
and.  Indeed,  appear  to  be  the  dominant  mobile 
entitles  In  the  conduction  process  occurring  In 
ordered  phases  of  compositions  near  p  “  0.50  at 
low  temperatures. 

II.  Computation  of  Ionic  Conductivity  Using 
the  MC  Method. 

The  conduction  plane  In  the  8 "-alumina  struc¬ 
ture  is  modeled  a*  a  two-dimensional  honeycomb 
lattice  In  which  all  sites  have  the  same  poten¬ 
tial  for  a  single  kind  of  mobile  cation,  and  a 
vacancy  mechanism  for  Ion  transport  Is  assumed 
(16).  We  choose  the  saddle-point  energy  for 
the  Na*  Ions  Co  be  U,  end  assume  that  two  Na* 
ions  on  nearest -neighbor  site*  loteract  repul¬ 
sively  with  an  energy  c.  No  Interaction  energy 
exists  between  a  Na*  Ion  and  a  NN  vacant  site 
or  between  two  neighboring  vacancies. 

The  Ionic  conductivity  of  our  model  for  Nag* 
alumina  1*  calculated  from  the  Nernst-Clnstelm 
equation! 

Ij  •  Sr2  •  <l> 

where  De  is  ths  charge  diffusion  coefficient,  C 
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He  have  performed  Langevln  dynamics  calculations  on  a  one-dimensional  model  for  frame¬ 
work  materials,  focusing  on  the  conductivity  of  interacting  mobile  Ions  In  the  pre¬ 
sence  of  Coulomble  craps.  We  observe  several  Interesting  phenomena,  arising  from 
competition  among  the  three  potentials  (lattice,  lon-lon,  ion-trap).  Findings  of 
Interest  Include:  L.  Generally,  Che  compensator  will  Indeed  trap  a  sufficient  number 
of  mobile  Ions  to  balance  Its  charge.  These  Ions  can  exchange,  but  they  are  quite 
well  localized  by  the  trap.  2.  For  Incommensurate  moblle-lon  concentrations,  the 
overall  DC  conductivity  Is  decreased.  3.  For  commensurate  concentrations,  the  DC 
conductivity  can  actually  Increase,  as  the  trapping  renders  the  remaining  moblle-lon 
array  effectively  Incommensurate. 


I.  Introduction 

Framework  solid  electrolytes  nearly  always 
require  disorder  (high  concentration  of  defects) 
on  one  sublattice  for  high  conductivity  to 
occur.  Ordinarily,  this  Is  achieved  by  attain¬ 
ing  a  formally  non-stolchlomctrlc  composition, 
such  as  (Na20)  1^.11*1203  for  8  alumina.  The 
concent  ration  of  mobile  ions  can  often  be  varied 
by  Introduction  of  formally  electrlcally-neutral 
components  (such  as  Sa20  above)  Into  the  struc¬ 
ture.  If  these  are  thought  of  as  fully  disso¬ 
ciated,  then  the  carriers  (Na+)  are  accompanied 
by  charge-compensating  counterions  (0"),  which 
remain  in  the  structure,  usually  at  rather 
clearly  well-defined  positions,  and  can  there¬ 
fore  effect  mobile  ion  motion  in  three  ways:  by 
distortion  of  the  framework  lattice,  by  blocking 
of  pathways,  and  by  Coulomb  trapping.  Theo¬ 
retical  studies  [l]  of  Ionic  motion  In  framework 
electrolytes  have  largely  Ignored  the  effect  of 
charge-compensating  Ions,  and  assumed  either  a 
simple  lattice-site  model  or  a  simple  periodic 
background  due  to  the  framework.  The  cole  of 
charge  compensators  has  been  studied  by  the  Oak 
Ridge  group,  who  present  evidence  for  motions  in 
the  spinel  block  due  to  Coulomble  forces  [21, 
and  more  specifically  by  Wolf  [31,  who  intro¬ 
duced  a  mechanistic  transport  model  In  which 
very  strong  trapping  by  compensators  Is  Invoked 
to  discuss  the  temperature  dependence  of  the 
Haven  ratio. 

More  general  considerations  were  presented  long 
ago  by  Mott  [41,  the  context  of  Che  mstal-lnsu- 
lator  transition.  He  argued  that  the  Coulomble 
crapping,  which  is  strong  for  low  concentrations 
of  charge  carriers,  will  become  screened  as  the 
carrier  concentration  increases,  snd  that,  when 
the  screening  length  Is  of  the  same  order  as  the 
lattice  constant,  the  formerly  trapped  charges 
become  free,  and  the  (exclconlu)  Insulator 
becomes  conductive.  Similar  ldeaa  seen  reason¬ 
able  for  ionic  materials,  though  the  screening 
effect  should  be  reduced  somewhat;  It  should  be 


more  reminiscent  of  Debye-fluckel  screening  In 
liquid  electrolytes  [5) 

Several  charge  coapensaclon  mechanisms  have  been 
postulated  for  noe-stolchlometrlc  Ionic  conduc¬ 
tors.  Three  possible  mechanisms  for  the  compen¬ 
sation  of  excess  cations  are:  1)  cations  in  the 
framework  are  replaced  with  cations  of  Lower 
oxidation  state.  For  example.  In  NASIC0N 
(Na3Zr2PSl20i2),  the  Sa*  Ion  concentration  can 
be  Increased  when  the  extra  positive  charge  is 
balanced  by  replacing  «tth  Sl+4  In  the 
framework;  2)  extra  cations  can  be  placed  in  the 
conduction  region  If  there  exists  a  correspond¬ 
ing  number  of  cation  defects  in  the  framework; 
and  3)  charge  neutrality  can  be  maintained  if 
the  excess  cations  are  balanced  by  interstitial 
anions. 

Trapping  sites  will  be  generated  by  these  three 
charge  compensation  mechanisms  If  the  deviations 
In  the  framework  structure  are  localized  near 
the  conduction  pathway.  Trapping  sites  might 
also  be  created  by  impurities  In  the  framework 
not  associated  with  charge  compensation. 

Our  objective,  in  the  present  paper,  is  to  study 
the  competing  roles  of  ionic  correlations,  back¬ 
ground  periodic  potentials  and  compensator  trap¬ 
ping  on  the  conductivity  of  a  one-dimensional 
model  for  solid  electrolytes.  The  method 
employed  will  be  stochastic  Langevln  dynamics 
I6J,  which  we  have  used  previously  [7,dJ  to 
examine  the  role  of  correlations  on  the  fre¬ 
quency-  and  temperature-dependent  conductivity 
of  such  mscerlsls.  We  find  chat  there  Is  an 
Interesting  trade-off  among  the  three  potentials 
In  the  problem,  snd  that  Coulomble  trapping  can, 
under  special  conditions,  actually  increase  the 
conductivity. 

The  computations  are  performed  for  a  one-dimen¬ 
sional  model  consisting  of  charged  ionic 
carriers  moving  in  a  sinusoidal  potential  snd 
subject  to  trapping  by  off-sxls  charge  compeu- 
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Electron-hole  pair  excitations  upon  atom  impact  on  a  metal  surface  are  studied  in  a  framework 
of  a  one-dimensional  independent-electron  model.  The  method  employed  treats  electron  dynamics 
quantum  mechanically  and  the  atom  motion  classically,  and  the  two  are  coupled  through  the 
time-dependent  self-consistent  field  (TDSCFl  approximation.  A  variational  method  is  used  to 
calculate  the  time  evolution  of  the  electronic  wave  packet.  Calculations  were  carried  out  for  the 
colliders.  He.  A r  and  H;  the  surface  parameters  were  chosen  to  model  Li.  Some  of  the  results 
obtained  are:  (1)  Electron  excitation  by  H  is  much  more  efficient  than  for  a  rare-gas  collider. 
Experimental  search  for  hole-pair  excitations  should  thus  be  best  pursued  with  H  as  a  collider.  (2) 
At  0  K  surface  temperature  <!£/£,  the  fraction  of  collision  energy  converted  to  hole-pair 
excitations,  decreases  as  the  collision  energy  increases  for  energies  up  to  - 1  eV.  At  collision 
energy  £  «  0.01  eV,  the  fraction  of  energy  transferred  is  -  0.2%  for  He  and  - 10%  for  H.  (3) 
Atom  trapping  due  to  energy  transfer  to  electrons  occurs  with  high  probability  (50-100%)  at 
sufficiently  low  collision  energies.  Ar  trapping  takes  place  at  energies  below  1  K  and  ft  trapping 
below  20  K.  (4)  The  calculations  show  a  pronounced  transition  from  atom  de-excitation  to  atom 
excitation  by  electron-hole  pairs  as  surface  temperature  increases.  (5)  Perturbation  theory  is  tested 
against  the  present  method.  It  breaks  down  mainly  for  trapping  and  for  temperature  effects. 
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Abstract 


A  Monte  Carlo  method  is  developed  for  simulation  of  mixed  ionic  conduc¬ 
tivity  in  8"-aluraina-type  materials.  The  conduction  plane  of  these  materials 
is  represented  by  a  lattice  gas  model  in  which  monovalent  and  divalent  cations 
carriers  diffuse  via  a  vacancy  mechanism  and  interact  through  a  nearest- 
neighbor  Coulombic  repulsion. 

By  comparing  experimental  data  for  pure  Na+  and  pure  Ba++  8"-aluminas  with 
simulation  results,  it  is  possible  to  estimate  the  near-neighbor  interaction 
energies  and  jump  barriers  Uj  for  both  kinds  of  ions.  On  the  basis  of 
these  estimations  the  total  ionic  conductivity  of  Na+-Ba++  3'*-alumina  is 
calculated  as  a  function  of  temperature  and  concentration  of  carriers.  As 
Ba4"4"  replaces  Na+,  the  conductivity  initially  increases  as  more  vacancies 
become  available.  For  very  high  temperatures,  thi3  increase  continues  until 
exchange  is  complete;  but  at  lower  temperatures,  the  conductivity  reaches  a 
peak  for  some  optimal  Ba++/Na+  composition,  and  then  drops  off  as  the  number 
of  Ba4"4",  and  hence  the  strength  of  ionic  correlation,  goes  up.  The  presence 
of  ordering  in  the  fully-exchanged  (all  Ba4-4")  case  manifests  Itself  in 
substantial  curvature  of  the  Arrhenius  plots  for  conductivity.  The  activation 
energy  for  conductivity  as  a  function  of  Ba4-4"  mole  fraction  (xBa++)  shows  a 
pronounced  rise  near  an  X  value  of  2/3,  in  agreement  with  recent  experimental 


observations 
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